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Abstract 

The electronic transport properties of a metallic carbon nanotube with the five-seven discli- 
nation pair characterized by a lattice distortion vector are investigated. The influence of the 
disclination dipole includes induced curvature and mixing of two sublattices. Both these factors 
are taken into account via a self-consistent perturbation approach. The conductance and the 
Fano factor are calculated within the transfer-matrix technique. 

PACS: 73.63.Fg, 72.80.Rj, 72.10.Fk 

INTRODUCTION 

Transport properties of variously shaped carbon nanostructures and graphene are of great prac- 
tical and theoretical interest. In particular, the conductivity of carbon nanotubes is currently 
the subject of wide investigations. Actually, there are many types of defects (such as vacancies 
and vacancy pairs [I], adatoms [2], structural defects, etc.) which influence the conductivity of 
carbon nanotubes. For instance, the electronic and transport properties of nanotubes containing 
vacancy pairs were found to be extremely sensitive to the sublattice positioning of the pair [3]. 
Moreover, it was observed that they substantially affect the metallic or semiconducting charac- 
ter of the tube. The structural (topological) defects in carbon nanotubes are mainly five-seven 
disclination pairs [3] and the Stone- Wales defects ("quadrupoles" , 5-8-5 or 5-7-7-5 defects |5|). 
The Stone- Wales structural defects were found to influence at least the local electronic proper- 
ties of the nanotube. It is also known that disclinaton topological defects (fivefolds) can convert 
planar graphene surface into the conical one with a marked modiflcation of the electronic states 
[7J. In the nanotubes, however, the presence of an isolated pentagonal ring is not allowed and 



the simplest topological defect is the 5-7 disclination dipole (DD). The simple model for the 
junction between two metallic tubes was investigated in [8]. 

In this paper, we study the electronic transport properties of a metallic nanotube containing 
closely spaced 5-7 disclination dipole (i.e., the junction composed of two semi-infinite metallic 
nanotubes). The paper is organized as follows. In Section [U the geometrical characteristics 
of the nanotube with the DD-induced curvature are discussed. The structure of a junction 
consisting of two nanotubes with different diameters is presented. In order to design the shape 
of the tube we introduce two parameters: the distortion vector b which influences both the 
radius and the chirality of the tube, and the phenomenological parameter characterizing the 
size of the curved region due to elastic relaxation. In Section [2l the effective Dirac equation is 
formulated, and a standard transfer matrix method is used to calculate the conductance and 
the Fano factor of the structure. The obtained results are discussed in Section [2l 

1 GEOMETRY 

A metallic or semiconducting character of the tube is governed by the translational vector 
To(n,m), where (n,m) are the numbers of steps along two unit cell vectors in the honeycomb 
lattice. The boundary conditions for the tube with a translational vector T transforms to the 
angular vector field a^p which dependis on chirality as (n + m) mod 3. This field, along with the 
angular momentum, subsequently generates the mass-type term in the effective Dirac equation. 
The value of the produced gap in the energy spectrum is determined as 2{j — a^)hVF/ Ro, 
where j is the angular momentum and Rq is the tube radius. When this term vanishes, the 
metallization of the nanotube occurs for j = (see [TO] for detail). 

As for the fivefold-sevenfold pair in the nanotube, it appears to be the source of translational- 
type holonomy [6j, therefore it should produce additional gauge vortex field. However, in this 
paper we restrict our consideration to a special case when the condition (n + m) mod 3 = is 
fulfilled at both sides. In other words, both tubes far from the junction region are suggested 
to be metallic. The gauge field produced by the dipole source should affect the chirality but 
the value (n + m) mod 3 remains untouched. Additionally, at low energies below the first gap 
value (|£'| << hVp/Ro, where Vp is the Fermi velocity and R is the tube radius) one can take 
into account for the estimation of first-order perturbation only the main conducting channel 
corresponding to the lowest angular momentum. 

As is known, the dislocation-type defect touches both the chirality and the radius of the 
carbon nanotube. In the general case, the tube with DD can be characterized by the distortion 
vector b: for the translational vector T one can find, that it has the value Tq, Tq = IttRq on 
the one side and Tq -|- 6 on the other. 

Trying to describe a shape of the nanotube, it is useful to perform a development of the 
junction region onto a 2D plane (see Fig. [T|). 

Generally, the DD with a distortion vector b changes both the length and the direction of the 
translational vector. As a result, on the right side the shape of the tube is also a cylinder with 
the radius \Tq + 6|/(27r), displaced by the angle between vectors Tq and Tq + b. In this paper, 
we consider b/To as a small parameter. As an example, this could be the vector 6(2, 1) which 
preserves the metallicity of the tube. Reverting to the shape of the tube we should roll up the 
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Figure 1: The schematic structure of the nanotube in the defect region. The fivefold is located at 
the point Q (Q'); and the sevenfold - at the point P (P')- Both the length and the direction of the 
translational vector Tq are changed in the presence of the 5-7 disclination dipole. 



development. This resembles the " cut-and-glue" procedure for a disclination in the nanocone 
[7J. This procedure determines an explicit relative positions for the cylinders (both the angle 
and the distance between the axis lines) as a function of b. For the exact shape of the structure 
in the junction region, let us construct a phenomenological approximation, which satisfies all 
the boundary conditions described above. 

Let us associate the development with the xy-plane and choose the axis along the tube on 
the left side as the x-axis. In this case, one has ro(0, ro,0) and r(^) = T{Tx,Ty,0). It is 
convenient to introduce the new frame with the vector n = {T x ez)/T {ei (i = x,y,z) are the 
orthogonal basis vectors) which follows the tube axis (see Fig. [2]). 

The orthogonal to n vectors are ri = — and T2 = ti x fi = T/T. To take into account 
the disclination dipole, which is situated at the top of the cylinder, one needs to describe the 
skewed conical surface between Q and P points. Actually, the DD has two effects: it changes 
the radius and shifts the surface of the tube in the direction normal to the tube axis. Since the 
radius R{S,) = T(^)/(27r) changes its value on AR, the ring will be shifted in the direction of T2 
by AR, so that we must add the shift T^Tji^i:). 

Finally, the surface is parametrized as following: 

T T T 
f = n^ + f2— + fi— cos(^ + f2— sin(^, (1) 
zvr zvr zvr 

where — oo < ^ < oo and < 99 < 27r are the normal and the transversal coordinates, respec- 
tively. 

The form of the tube with 5-7 rings is supposed to be 
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Figure 2: A schematic picture of the nanotube surface. The fivefold is located at the point Q, and 
the sevenfold - at the point P. The thick tube axis BB' coincides with the x-axis. The thin tube axis 
CC is shifted due to the disclination dipole QP. 



Figure 3: A molecular-dynamics simulation for the shape of the (12,0)-(10,1) junction [9]. Bottom: 
in the junction region the shifting of the tube axis can be observed. 
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where we assume 



n6=ro + ^(tanh^ + l) 



(3) 



to be the tube translation vector, depending on the coordinate ^, the distortion vector b and 
an effective parameter which corresponds to the half-width of the surface region curved by 
the DD. Since the surface of the nanotube has high but finite Young modulus, one can expect 
that .^0 is of the same order as Rq. Let us calculate the geometrical properties of the tube in 
the first order in b. 

The metric tensor is found to be 



5ee = l+e 
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g^^ = 1- 7^^, g^f 
The metrical connection 
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is found to be 



r|e = %€^/^o-5^7ee/(2^o); 

= = ^a^^/(2i?g), 

and = 0. Tetradic coefficients e*^, which are determined by the equation 



are chosen to be 
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The spin connection coefficients 

where := 3^ + is a covariant derivative, reads 
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2 CONDUCTIVITY AND THE FANO FACTOR 

From now on we set Ep = 0, Vp = h = 1. The Dirac equation on the curved surface is defined 



as 



-ia^e>'{V^,-iA^)lP = ElP, (12) 

with (Tj being the Paufi matrices {i = 1,2), -0 is a 2-component spinor wavefunction, and 
(// = ^, (f) being the DD field. The derivative includes a spin connection term Jl^ and it is 
written as 

V^ = d^ + n^, = -uj^''[aa, (Jb]. 
In the linear in b approximation the Dirac equation takes the form 

- iaie^^ (Vg - iA^)i> - ia2{e^^di: + e^^^iV^ - iA^))i> = E-iP, (13) 
with the covariant derivatives written as 



167r2i?o 

n/f 

Xl^ 

87r2i?o ' ^ 47r2i?2 J 2 

Substituting p4p into (\l'd^ we find both the undisturbed Dirac operator T>o 



= ^« + [-^TsTT + g .^2p2 COS v^. (14) 



Po = -idiS^ - i^d^, (15) 
-Ko 

and the perturbation operator V 

V = iai'^d^ - iaifl^ - aiA^ + ^'^2^9^ - - + ^'^2^^d^- (16) 

It should be mentioned that the perturbation (|16|) includes the last term, which is not 
localized (differs from zero value at ^ — )• 00). Therefore an analysis of (jlSp and ()16p differs from 
a typical scattering problem. To simplify the problem, let us restrict our consideration to the 
lowest term with the zero angular momentum, that is only a plain wave ip{^,ip) = (1, zbl)"^e*'^^ 
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is considered. In this case in (jlSp and (|16p only the zero-momentum terms should be taken into 
account in the Fourier series for the angle-dependent part of the wavefunction. 

As a result, instead of (jlSp and (jl6p there should apper an equation for the zero-momentum 
part of the wavefunction with the Dirac hamiltonian Dq = —iaid^ and the perturbation 

Here <><^ denotes the averaging over if and (6xTo)^ = bxTQy—Toxby. Now the perturbation (fT7|) 
is localized near ^ = with a half- width ^q. In order to remove the differentiation in (jl7p one 
can multiply the operator Dq + V by 

A ^ ibx To), \ 



87r2i?oeocosh2(e/eo); ■ 

In this case, the additional potential 

IE{bxfo)J{87r^Ro^ocosh\^/^o)) 

appears, where / is a unit matrix. 

The field A describes the DD-induced valley mixing. It may be written as Aj = ^ej^d^G 
where j, A; = ^, is an effective "charge" of each disclination, G is a dipole potential, which 
satisfies the equation 

27]- 

= d'^ + d'^/ Rq, Ejk is a unit antisymmetric tensor, 5[x) and 6'{x) are the Dirac delta function 

and it's derivative, and d is a vector connecting five- and sevenfold (it has the same order as h 
Explicitly, one can find G in a form 

G{i,^) = -'^Y.id^signiO + ^sign(n)]e-'^e-l<l/^«. (18) 

RO n ^ 

For the zero momentum n = 0, the field A has the form 

= 0, = ^^5(0, (19) 

where (5(^) is a Dirac delta function. Finally, the perturbation operator looks as following: 

r iby<fo),E b-fo 47rV6;g 

87r2i?oeocosh2(e/eo) 'l6^2^g^ocosh2(C/eo) ' Rl ^ ^ 

Notice that the last term in (|2Up corresponds to the .^-dependent mass term cT2m(.^) in the 
Dirac equation. Let us calculate the transfer matrix and the conductance according to the 



7 



method described in [12] and [13J. To transform the initial Dirac equation to the diagonal form, 
the unitary rotation is performed: V ~^ ^V') 1^ ~^ CVC\ ai — )• CaiC^ , C = {ai + as)/V2. For 
zero angular momentum, the equation for the transfer matrix reads 

r(f ,-^) = T^'\^,-^)-^ j\ dxT^^\^,x)a-Mx)Tix,-^), (21) 

2 

where T^'^\x2,xi) = exp{ia3E{x2 — xi)) and L is the length of the tube. The perturbation ([20]) 
takes the form 

(73 V = -0-3— —— . o ,^ ,^ , + J - 



'87r2i?oeocosh2(C/eo) e + ^^Rl 

-'l6,=fl§{„cosh^e/£„)-""-Rr**^'- 

We assume that L is much larger than ^o- In the linear in h approximation one can replace T 
by 7''-''^ in the integral in ()2ip . In our consideration, the direct intervalley scattering is not taken 
into account. Thus one can consider two inequivalent i^-points (valleys) as two different channels 
coupled by the non-Abelian field. Therefore, one should perform a replacement /i — )• with 
jjL = 1/A and T2 being the Pauli matrix in the channel space representing the mixing of sublattices 
due to pentagonal and hexagonal defects (see [7j for detail). By substituting ([22]) into (f2T]l and 
taking the limits of integration to be infinity, one obtains the transfer matrix in the form 

7 = exp(z(T3£;L) + ^[ ^^2ji2 + ^^3 4^2fig — J exp(zo-3^L) + fid^aiT2^. (23) 

The scattering matrix (t^)~^ is defined by the upper-left element of the transfer matrix, so that 
the last term in ()23p is neglected. The conductance is defined as G = (4e^//i)rr[t^t], where the 
trace is taken over the channel space. It takes the form 

where e = ERq is a dimensionless energy (in the absolute units of hVp/Ro). Correspondingly, 
the Fano factor F =1- Tr[{thf]/Tr[th] reads 

One can see that the mass-type term (the DD term) does not influence both the conductance 
and the Fano factor since only linear terms in b are taken into account (the mass is expected to 
be the second-order term). Depending on the angle between Tq and Tq + b, two types of energy- 
dependent behaviour appear. For the angle less than 7r/2, both the scalar and "vector" products 
in (|24p and (|25p are of the same sign. In this case, the linear in energy term in (|24p leads to a 
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Figure 4: Conductance G, in the units of hj (4e^), depending on the energy E in the units of hVp/Ro, 
through the metalhc junction (9,0) - (11,1) tubes (black hne) and (9,0) - (10,2) (gray hne). 

decrease of the conductance with increasing energy (and vice versa), see FigH] (black line). This 
case corresponds to slightly curved tubular structures. For the second type of behaviour, when 
the angle is greater than 7r/2 the increasing of conductance is expected with energy increasing. 
A special case is presented on FigU] (gray line), with Tq and b orthogonal, where a transmission 
coefficient is equal to unity at the Fermi energy e = 0. However, in this case the higher-order 
perturbation terms should be included. 

CONCLUSION 

In this paper, we have investigated the electronic transport properties of metallic nanotubes 
with the disclination dipole defect. The case of disclination dipoles preserving the chirality 
(metallic character of the tube) was considered. This assumption simplifies the problem under 
consideration and allows us to construct a self-consistent expansion scheme for the effective Dirac 
equation. The perturbation includes both the curvature term, which appear as the Gaussian-like 
smooth localized potential, and the DD-induced term. The standard transfer-matrix approach 
is used to calculate the conductance and the Fano factor of the structure. 

As the marked effect of the curvature, the energy-dependent term in both the conductance 
and the Fano factor appears. It should be stressed that the tube length L and the effective 
half-width do not enter the final results. This differs from the case of ballistic transport in 
graphene with various types of disorder studied in [T3], where both the conductance and the 
Fano factor were found to depend on the sample length (or, more precisely, on EL). 

Notice that in [14], ^ a similar model for two metallic nanotubes of different radiuses con- 
nected by the conical interface was constructed. The transmission probability was found to 
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decrease when the difference of tube radiuses (which is determined by by in our paper) in- 
creases, which is in agreement with the 1 — ab'^ behaviour in (|24p . The origin of this behaviour 
within our model is purely geometrical, while in |14] the boundary conditions for the wave- 
function incorporates both the "geometrical" factor and the mixing of the K-points within the 
junction. As for the energy dependence, in [8] the conductance was found to decrease with 
energy increasing for (17,17)-(18,21) junction, and subsequently the increasing of conductance 
was found for (23,8)-(16,22) junction. This fact is in general agreement with ()24p . where the 
character of energy behaviour depends on the angle between translational (chiral) vector and 
the distortion vector (the difference of chiral vectors). At higher energies close to band edges, 
a more complex behaviour of conductance was predicted, which is absent in our model. The 
source of this difference is twofold. First, instead of the boundary conditions on the tube-cone 
interface in |14) and^ we consider the model with smooth shape of the structure, so that the 
resonant behaviour of wavefunction in the junction region does not appear in our case. Second, 
for simplicity terms with higher momenta were not included in our model and, as a result, lo- 
calized and resonant states are absent even at low energies. A self-consistent inclusion of higher 
momenta in our approach is an important open problem. 

The disclination dipole is taken into account via the non-Abelian field, which leads to the 
mixing of valleys in the low-energy behaviour. The infiuence of this field was found to be 
negligible within the perturbational approach compared to the infiuence of curvature, in contrast 
with the case of ordinary vortex topological field sources [15]. A dominating process of valley- 
mixing was taken into account (as it is estimated for the scattering in carbon nanocones [7]) 
instead of the real intervalley scattering, which is expected to be omitted. This results in the 
correction of the conductance and the Fano factor in the highest order in b to have the same 
form as for the single-channel approximation. One should also note, that within our model the 
chirality-dependence of the potential was neglected, because both the left and right ends of the 
tube with the disclination dipole are of the metallic-type (despite the fact that the radius and 
exact chirality indices are different). In the general case, one should expect the infiuence of the 
defect on transport properties to be sensitive on the tube chiralities, as it was observed in [3]. 
Another important question is the finiteness of the free electron path in the real tube. As it 
is easy to see, the phenomenological localization length of the curved area does not infiuence 
the transport properties within the constructed approximation. Due to this fact, an effective 
long-range infiuence of topological defects could possibly lead to the non-trivial interplay with 
the mean free path length in the tubes with disclination dipoles. 
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